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ABSTRACT 

We have analyzed ~600 ks of Rossi X-ray Timing Explorer data of the neutron star low-mass X-ray 
binary and Z source GX 17+2. A study was performed of the properties of the noise components and 
quasi-periodic oscillations (QPOs) as a function of the broad-band spectral properties, with the main goal 
to study the relation between the frequencies of the horizontal branch (HBO) and upper kHz QPOs. It 
was found that when the upper kHz QPO frequency is below 1030 Hz these frequencies correlate, whereas 
above 1030 Hz they anti-correlate. GX 17+2 is the first source in which this is observed. We also found 
that the frequency difference of the high frequency QPOs was not constant and that the quality factors 
(Q values) of the HBO, its second harmonic, and the kHz QPOs are similar, and vary almost hand in 
hand by a factor of more than three. Observations of the normal branch oscillations during two type I 
X-ray bursts showed that their absolute amplitude decreased as the flux from the neutron star became 
stronger. We discuss these and other findings in terms of models that have been proposed for these 
phenomena. We also compare the behavior of GX 17+2 and other Z sources with that of black hole 
sources and consider the possibility that the mass accretion rate might not be the driving force behind 
all spectral and variability changes. 

Subject headings: accretion, accretion disks - stars: individual (GX 17+2) - X-rays: stars 



1. INTRODUCTION 

Based on their broad-band spectral and variability prop- 
erties, six of the persistently bright neutron star low- 
mass X-ray binaries (LMXBs) were classified as Z sources 
(Hasinger & van der Klis 1989), after the Z-like tracks 
they trace out in X-ray color-color (CDs) and hardness- 
intensity diagrams (HIDs). These sources are GX 17+2, 
Cyg X-2, GX 5-1, GX 340+0, Sco X-1, and GX 349+2. 
The Z-like tracks consist of three branches, which, from 
top to bottom, are referred to as the horizontal branch 
(HB), the normal branch (NB), and the flaring branch 
(FB). It is generally believed that the parameter that de- 
termines the position along the Z track is the mass accre- 
tion rate, increasing from the HB to the FB. Note that 
the definition of the mass accretion rate in this respect is 
rather vague - it is often defined as the mass accretion rate 
through the inner disk and onto the neutron star surface. 
In addition to spectral changes along the Z-track, some of 
the Z sources show long term changes in the shape and 
position of the Z-track in the CD and HID. These secu- 
lar changes, as they are referred to, are clearly observed 
in Cyg X-2 (Kuulkers et al. 1996; Wijnands et al. 1997b), 
GX 5-1 (Kuulkers et al. 1994), and GX 340+0 (Kuulkers 



& van der Klis 1996) (the Cyg-like sources; the other Z 
sources are referred to as the Sco- like sources), and more 
recently also in GX 17+2 (Wijnands et al. 1997a). It has 
been suggested that they are related to the relatively high 
inclination at which these sources are seen (Kuulkers et al. 
1994; Kuulkers & van der Klis 1995), or to different mag- 
netic field strengths of the neutron stars (Psaltis et al. 
1995). 

The power spectra of the Z sources show several types of 
quasi-periodic oscillations (QPOs) and noise components 
(see van der Klis 1995a, for a review). It was found that 
their presence and properties are very well correlated with 
the position of the source along the Z track (Hasinger & 
van der Klis 1989), even when the Z-tracks show secu- 
lar changes (e.g. Kuulkers et al. 1994). Three types of 
low frequency (<100 Hz) QPOs are seen in the Z sources: 
the horizontal branch (HBOs), normal branch (NBOs) and 
flaring branch QPOs (FBOs). Their names derive from 
the branches on which they were originally found. The 
HBO is found on the HB and NB with a frequency (15- 
60 Hz) that gradually increases along the HB towards the 
NB. When the sources move from the HB onto the NB 
the frequency increase flattens off. In GX 17+2 and Cyg 
X-2 it was found that when the source passes a certain 
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point on the NB, the HBO frequency starts to decrease 
(Wijnands ct aL 1996, 1997b). The NBO and FBO are 
most likely the same phenomenon. They are found on the 
NB and FB (near the NB/FB vertex) but not on the HB. 
On the NB the QPO has a frequency of ^5 7 Hz, which 
rapidly increases to ~20 Hz when the source moves across 
the NB/FB vertex (Priedhorsky et al. 1986; Dieters & van 
der Khs 2000). In recent years two high frequency (or 
kHz) QPO were found in the Z sources (van der Klis et al. 
1996; Wijnands et al. 1997a, 1998a,b; .Jonkcr ct al. 1998; 
Zhang et al. 1998, see van der Klis 2000, for a review). 
They are often observed simultaneously, with a frequency 
difference of ^^300 Hz, and have frequencies between 215 
Hz and 1130 Hz, which increase from the HB to the NB. 
In Sco X-1 the frequency difference was found to decrease 
with increasing QPO frequency (van der Klis et al. 1997b). 
The frequency difference in the other Z sources is both con- 
sistent with the behavior seen in Sco X-1 and with being 
constant (Wijnands et al. 1997a; Jonkcr ct al. 1998; Psaltis 
et al. 1998). Three types of noise are seen in the Z sources. 
They are the very low frequency noise (VLFN), the low fre- 
quency noise (LFN), and the high frequency noise (HFN). 
The VLFN and HFN arc found on all branches, whereas 
the LFN is only observed on the HB and NB. The VLFN, 
which is found at frequencies below 1 Hz, can be described 
by a power law. The HFN and LFN are both band lim- 
ited components, with cutoff frequencies of, respectively, 
10-100 Hz, and 2 10 Hz. 

Many competing models have been proposed for the ori- 
gin of the QPOs and noise components. It is beyond the 
scope of this introduction to mention these models in de- 
tail - most of them will be discussed in Section 4. 

In this paper we present a study, based on data acquired 
with the Rossi X-ray Timing Explorer {RXTE), of the cor- 
related spectral and variability properties of the Z source 
GX 17+2. It is a continuation of the work by Wijnands 
et al. (1996, 1997a). The current paper constitutes the 
first report on the very large observing campaign of 1999, 
which more than doubled the total coverage of the source. 
This campaign was undertaken with the express purpose 
of investigating if a non-monotonic relation exists between 
the frequency of the kHz QPOs and the HBO in GX 17-^2. 
Section 2 deals with the observations and analysis. The 
spectral results are presented in Section 3.1 and the re- 
sults for each power spectral component in Sections 3.2 
and 3.3. A number of qualitatively new results is found 
in our greatly expanded data set. In particular, we find 
that when the upper kHz QPO frequency is below 1030 
Hz, it correlates with the HBO frequency, whereas above 
1030 Hz they anti-correlate. We also find that kHz QPO 
frequency difference is not constant and that the Q values 
of the HBO and kHz QPOs are probably not explained 
by life time broadening. These and other results are com- 
pared to observations of other Z sources and discussed in 
terms of current models in Section 4. Finally, we argue 
that changes in the spectral and variability properties of 
the Z sources may not be not driven by changes in the 
mass accretion rate. 

2. OBSERVATIONS AND ANALYSIS 

The data used for the analysis in this paper were all 
obtained with the Proportional Counter Array (PCA, Ja- 



hoda et al. 1996) on board RXTE (Bradt ct al. 1993). 
The PCA consists of five xenon-filled proportional counter 
units (PCUs), each with an effective area of ~1250 cm^ (at 
10 keV). Although the five PCUs are in principle identical 
they all have a slightly different energy response. These 
responses change continuously due to slow processes such 
as gas leakage and the aging of the electrodes. In addition, 
the high voltage settings of the instruments are occasion- 
ally altered (gain changes), resulting in rather more drastic 
changes in the detector response. These changes have, at 
the time of writing, been applied three times during the 
life time of RXTE thereby defining four gain epochs. Oc- 
casionally one or more PCUs are not operational. They 
can be switched off by an internal safety mechanism, or 
by the ground control crew, for reasons of detector preser- 
vation. Therefore, the number of active detectors varies 
between the observations. 

All our RXTE/PCA observations of GX17+2 were done 
between 1997 February 2 and 2000 March 31. A log of the 
observations is given in Table 1. We do not include the 
observations done in February 1996, which were used by 
Wijnands et al. (1996). The reasons for this are the lim- 
ited time resolution in the energy range of interest and 
difficulties with scaling the Sz (see below) due to an in- 
complete Z track. Data taken during satellite slews and 
Earth occultations were removed. The nine type I X-ray 
bursts that were observed are the subject of a separate 
article by Kuulkcrs ct al. (2001); two of those were also 
studied in this paper. The total amount of good data that 
remained was ~600 ks. 

Data were collected in several modes with different time 
and spectral resolutions. Two of these modes, 'Standard 
1' and 'Standard 2', were always operational. In addition 
to these two modes, other modes were active that varied 
between the observations. The properties of all modes are 
given in Table 2. 

The Standard 2 data were used to perform a broad-band 
spectral analysis. The data were background subtracted, 
but no dead time corrections were applied; these were in 
the order of 2-5%. 

For each 16 s data segment (i.e. the intrinsic resolu- 
tion of the Standard 2 mode) we defined two colors, which 
are the ratios of count rates in two different energy bands, 
and an intensity, which is simply the count rate in one 
energy band. The energy bands used for the colors (soft 
and hard color) and intensity, are given in Table 3. The 
lower energy boundaries for the soft color and intensity 
were chosen relatively high in order to avoid the lowest 
and least reliable energy channels. By plotting the two 
colors against each other a color-color diagram (CD) was 
produced. A hardness-intensity diagram (HID) was pro- 
duced by plotting the hard color against the intensity. To 
produce the CDs and HIDs we only used data obtained 
with those PCUs that (for each gain epoch) were always 
on. For gain epoch 3 these were PCUs 0, 1 and 2, and for 
gain epoch 4 PCUs and 2. Due to the different numbers 
of detectors and the differences in the detector settings we 
decided to produce the CDs and HIDs separately for the 
two gain epochs. However, in choosing the energy chan- 
nels we tried to take channels whose energy boundaries 
were as close as possible. 

Due to the aging processes mentioned earlier, observa- 



3 



tions with the same source spectrum that are made more 
than a few weeks apart end up at a different location in 
the CD. To correct for this effect we analyzed a number 
of RXTE/PCA observations of the Crab pulsar/nebula 
(which is assumed to have a constant spectrum, sec also 
Kuulkers et al. 1994) that were taken around the time of 
our GX 17+2 observations. For all Crab observations we 
produced the colors in the same energy bands as we used 
for GX 17+2. We found that the observed colors of the 
Crab indeed changed. For each observation wc calculated 
multiplicative scaling factors, for both colors, with respect 
to those of the first Crab observation. We used these fac- 
tors to scale the colors of the GX 17+2 observations back 
to those of the first GX 17+2 observation. This procedure 
could only be applied to the epoch 3 observations, where, 
as expected, we found it to result in a narrower track in the 
CD. No corrections were applied for the intensity. Since 
no Crab observations were available for the March 2000 
observations, no corrections could be applied to epoch 4 
data. The shifts in that data set appeared to be small 
anyway. 

Our power spectral analysis was based on selecting ob- 
servations as a function of the position along the Z track 
in the CD. We used a method that is based on the 'rank 
number' and 'Sz' parameterization methods introduced by 
Hasinger et al. (1990) and Hertz et al. (1992), that has 
been gradually refined in similar studies (Kuulkers et al. 
1994; Kuulkers & van dcr Klis 1996; Kmilkcrs et al. 1997; 
Wijnands et al. 1997b; Dieters & van der Klis 2000). In 
this method in its current form, all points in the CD are 
projected onto a bicubic spline (see Press et al. 1992) 
whose normal points are placed by hand in the middle 
of the Z track (see Figure 1). The HB/NB and NB/FB 
vertices are given the values Sz=l and Sz=2, respectively. 
The rest of the Z track is scaled according to the length of 
the NB. Problems arise when sharp curves are present in 
the track, in our study most notably around the NB/FB 
vertex. Due to scatter, points which have FB properties 
end up on the NB and vice versa. This is not a problem 
that is only intrinsic to the Sz parameterization; it is a 
limitation that applies to all selection methods based on 
colors. No observable parameter, apart from the power 
spectra, has been identified that could be used to distin- 
guish NB and FB observations around the vertex better 
than the X-ray colors. 

In order to improve on the results of Wijnands et al. 
(1997a) we wished to combine the epoch 3 and epoch 4 
data sets. Unfortunately, the tracks traced out in the 
CDs of epoch 3 and 4 are not the same and two differ- 
ent splines had to be used for the Sz parameterization. 
Since the normal points for these splines were drawn by 
hand, the Sz scales for the epoch 3 and epoch 4 CDs 
were unlikely to be exactly the same. Therefore, we first 
transformed the Sz scale of epoch 4 to that of epoch 3. 
To accomplish this, we measured the frequency of either 
the HBO or NBO at several places along the Z track of 
epoch 3, and determined the Sz interval corresponding 
to the same frequency in the epoch 4 data. We found 

that Sz,epoch4 = (0.06 ± 0.04) + (1.005 ± 0.018)^,,epoch3, 

showing that, although the scales are the same (within the 
errors), a small shift is present. We subsequently scaled 
the Sz values of epoch 4 using the above expression. The 



above scaling method assumes that the HBO and NBO 
frequencies are strongly related to Sz ■ Previous studies of 
Z sources (see e.g. Kuulkers et al. 1997; Wijnands et al. 
1997a; Jonker et al. 2000b; Dieters & van der Klis 2000), as 
well as the fact that only little scatter was present around 
the linear relation seem to confirm this. 

Power spectra were created from the data in modes with 
high time resolution (< 2~^^ s; see Table 2) using standard 
Fast Fourier Transform techniques (see van der Klis 1989, 
and references therein). The data were not background 
subtracted and no dead time corrections were applied prior 
to the Fourier transformations. We made power spectra 
in several energy bands, with several frequency resolutions 
and Nyquist frequencies. We finally settled on 0.0625-4096 
Hz power spectra in the 5.1-60 keV (epoch 3) and 5.8-60 
keV (epoch 4) bands, since the QPOs were most signifi- 
cantly detected in those bands (note that the source does 
actually not contribute significantly to the power spec- 
tra above '^40 keV, however, our data modes did not al- 
low us to discard data above this energy). This choice of 
frequency range means that the properties of the VLFN, 
which dominates the power spectrum below 0.1 Hz, could 
not always be measured satisfactorily, but it does allow 
us to follow the power spectral evolution on time scales 
down to 16 s. Note that the power spectra of epoch 4 
were produced in a slightly higher energy band than those 
of epoch 3. Due to the gain changes and the limited choice 
of energy channels this was the best possible match. 

The power spectra were selected on the basis of Sz, as 
determined from the CD. The Sz selections usually had 
a width of 0.1 and did not overlap - no power spectrum 
was represented in more than one S'z-selcction. Differ- 
ent widths were used in cases where the power spcctrimi 
changed rapidly as a function of Sz (narrower selections) 
or when the powers were weak (wider selections). All the 
power spectra in a selection were averaged, and the re- 
sulting power spectrum was rms normalized according to 
a procedure described in van der Klis (1995b). 

The properties of the power spectra were quantified 
by fitting functional forms to them. The low frequency 
(0.0625 256 Hz) and high frequency (100 4096 Hz) parts 
of the power spectrum were fitted separately. The high 
frequency part was fitted with one or two Lorentzians 
(for the kHz QPO(s)) and with the function P(i') = 
Pi + P2Cos{2ttv/vm) + Pzcos{A'Ki' /vn) for the dead time 
modified Poisson level (Zhang 1995; Zhang et al. 1995); no 
separate term was used for the contribution of the Very 
Large Events (VLE) count rate since it was absorbed by 
this fit function. At low frequencies the noise that was 
present could not be fitted consistently, and varying com- 
binations of functionals had to be used: 

• S'^<0.0: A Lorentzian (LFN+VLFN) 

• 5^=0.0-0.1: A cut-off power law (LFN+VLFN) 

• 5^=0.1-1.4: A power law (VLFN) and a cut-off 
power law (LFN) 

• 5z=1.4-1.6: A power law (VLFN) and a Lorentzian 

(LFN+NBO) 

• 5z=l. 6-5.0: A power law (VLFN) and a Lorentzian 
(LFN) 
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• S'^>5.0: A power law (VLFN) 

The expression for a power law is P(j^) cx: i^^", that 
for a cut-off power law is P{i^) oc z/-«e~''/'^<="* (where Vcut 
is the cut-off frequency), and for a Lorentzian P{iy) oc 
l/[{u - + {FWHM/2)'^] (where is the centroid fre- 
quency and FWHM is the full- width-at-half- maximum). 
The HBO, its second harmonic and the NBO/FBO were 
each fitted with a Lorentzian. Below the HBO, at about 
half its frequency, a broad bump was present that was also 
fitted with a Lorentzian, whose frequency was sometimes 
fixed to zero. At low frequencies the dead time modified 
Poisson level was fitted with a constant. An example of a 
fit to a power spectrum (52=0.4-0.5), displaying the con- 
tributing of the several components, is shown in Figure 
4. 

Note that changes in the fit function, such as using a 
cut-off power law instead of a Lorentzian or adding an ex- 
tra component, may lead to minor changes in the values 
of other parameters. Errors on the fit parameters were 
determined using Ay^ = 1. Upper limits were determined 
by fixing some or all of the parameters of a component, 
except the rms amplitude, to values similar to those ob- 
tained in the closest Sz selection where it was found to be 
significant, leaving all other fit parameters free, and using 
Ax^ = 2.71 (95% confidence). 

A study of the energy dependence and time lag proper- 
ties of the QPOs and noise components will be presented 
elsewhere. 

3. RESULTS 

3.1. Broad-hand spectral behavior 

The CDs and HIDs of both epochs are shown in Figure 
1. In both CDs the HB/NB vertex is not well defined, 
since the HB is almost a continuation of the NB. This is 
mainly due to the relatively high energies that we chose 
for our soft color (sec Section 2). At lower energies the 
turn-over is much clearer. In Figure 2 we show the count 
rate in several energy bands as a function of S^, for epoch 
4. It shows that the HB/NB vertex in the HID is entirely 
due to the count rates at low energies. At high energies 
the HB is a perfect continuation of the NB, and no vertex 
is present. 

In the HID of epoch 3 the Z track appears to be seg- 
mented, which was already noted by Wijnands et al. 
(1997a). Although no corrections for the slow detector 
aging processes were applied to the intensity in this HID, 
we note that a quick look showed that they would only 
have made the segmentation more apparent. The shifts in 
the HID might be due to the secular motion that has been 
observed in other Z sources. The shifts, of up to ^^5% in 
intensity, do not show up in the CD. This is because colors 
are ratios of intensities and are therefore not very sensi- 
tive to overall intensity changes; exactly for this reason 
CDs are preferable over HIDs for our purpose. Moreover, 
the width of the tracks in the CD is about 5%, so changes 
smaller than this are hard to observe. 

The bottom panel of Figure 3 shows the distribution of 
the time spent by the source in each part of the Z track. 
The source spent 28% of the time on the HB {Sz<l), 44.2% 
on the NB {1<S:,<2), and 27.8% on the FB {Sz>2). The 
average speed along the Z track as a function of Sz{{Vz)) 



is shown in the top panel of Figure 3. The speed at a given 
SS) is defined as VS) = \Sz{i + 1) - Sz{i - 1)1/32 (see 
also Wijnands et al. 1997b), where i is used to number the 
points in order of time. As expected the (Vz) increases 
considerably when the source enters the FB, but also at 
the top of the HB [S^ < 0.5). Combined with the small 
amount of time spent in the upper HB we can conclude 
that the source reaches to Sz values this low only in brief, 
quick dashes. 

For a more detailed analysis of the spectral changes as 
a function of the position along the Z track in GX 17-1-2 
we refer to O' Brien et al. (in prep.). 

3.2. Power spectra 
3.2.1. Low Frequency QPOs and Noise Components 

Figure 5 shows the low frequency part (0.0625-256 Hz) 
of the power spectrum for nine different 5*2 selections. The 
power spectra shown in Figure 5 are selected from epoch 
3 and epoch 4 and are therefore a combination of 5.1-60 
keV and 5.8-60 keV data. The percentage of epoch 3 and 
epoch 4 data varies between the S'z-selections. The contri- 
bution of epoch 3 data is highest at the extremes of the Sz 
range, reaching 100%, and gradually decreases from both 
ends to ~10% around ^2=1. 

Several components can be seen, which change in 
strength and shape as a function of ^2 . The different com- 
ponents are identified in Figure 5 (see also Figure 4) and 
the ^2 ranges in which they were detected are given in 
Table 4. The component identified as the sub-harmonic of 
the HBO is rather broad and does not have the appear- 
ance of a QPO. However, based on the frequency ratios 
(see below) it is referred to as the HBO sub-harmonic. 

Some difficulties were experienced while fitting the low 
frequency part of the power spectrum between ^2=1. 4 and 
5*2=1.6, where two or more components with similar fre- 
quencies were simultaneously present. Above 5*2=1.4 the 
NBO appeared, as a broad feature on top of the LFN. We 
were not able to distinguish the two components and de- 
cided to fit them together with a single Lorentzian. The 
fit values are not used in figures and tables, since they do 
not represent any of the individual components. Above 
^2=1. 6, the NBO and LFN could be distinguished more 
easily. The fit function used for the LFN, which was un- 
derlying the NBO, was changed to a Lorentzian to be more 
consistent with fits at higher Sz (fits with a cut-off power 
law gave equally good Xred ^z)- 

In the following sections the results for each of the com- 
ponents will be presented. 

3.2.2. HBO 

The HBO was detected between 5*2=— 0.6 and 52=2.1 
and its second harmonic between ^2=— 0.6 and 5*2=1.0. 
The second and third columns of Figure 6 show their 
rms amplitudes, FWHM, and frequencies as a function 
of 5*2 (see also Table 5). The frequency of the HBO in- 
creased from 21.3 Hz at 5'2=-0.43 to 60.3 Hz at 5*2=1.45 
and then decreased to 48.2 Hz at 5z=2.04. The second 
harmonic of the HBO had a frequency that was on aver- 
age 1.941±0.007 times that of the first harmonic, which is 
significantly different from the expected value of 2. An ex- 
planation for this discrepancy is proposed in Section 4.1.2. 
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Between Sz=— 0.6 and Sz — ^-O both the HBO and its sec- 
ond harmonic decreased in strength, with the decrease of 
the second harmonic occurring somewhat faster than that 
of the HBO. In that same Sz interval the FWHM of the 
HBO and its second harmonic were fairly constant (show- 
ing a slight increase), although considerable scatter was 
present around the average values (which were, respec- 
tively, 7.9±0.3 Hz and 14.4±0.6 Hz). The Q values of 
the HBO and its harmonic are shown in Figure 7. In 
the Sz range whore both were detected their Q values 
were consistent with each other. When the second har- 
monic was not significantly detected anymore {Sz>1.0) 
the rms amplitude and the FWHM of the HBO both 
decreased from, respectively, 2.45±0.08% and 9.7±7 Hz 
{Sz=0.5 1.0) to 1.48±0.06% and 5.1±0.4 Hz (S'^^1.0- 
1.5). Between Sz=l-0 and 52=1.5 the relation between 
the Sz and the HBO frequency started to flatten. Above 
S'z^l.5 the frequency of the HBO clearly dropped, ini- 
tially quite abruptly and later on more smoothly. This 
frequency drop coincided with an increase in the FWHM 
to 15.7±1.5 (52=1.5-2.1); the rms amplitude showed a 
small increase to 2.1±0.1% (5*2=1.76), followed by a de- 
crease to 1.2±0.2% (52=2.04). 

Underlying the HBO and its second harmonic, but with 
a central frequency lower than that of the HBO, we found 
a broad feature that was also fitted with a Lorentzian. It 
was significantly detected between 52=— 0.6 and 52=1.3. 
The properties of the broad feature arc shown in the first 
column of Figure 6 (see also Table 5). Below 52=0 the 
frequency of the Lorentzian was fixed to zero. Between 
52=0.5 and 52 = 1.5 the frequency of the broad feature 
is on average 0.539±0.015 times that of the HBO. This 
could mean that the broad feature is the sub-harmonic 
of the HBO, certainly when one takes into account that 
the frequency of a broad feature is rather sensitive to the 
shape of the continuum. Whereas the rms amplitude and 
the frequency of the sub-harmonic both change strongly 
with 52 the FWHM remains more or less constant, with 
an average value of 29.5±1.5 Hz. 

3.2.3. NBO and FBO 

The NBO and FBO were detected between 52 = 1.6 and 
52=2.7. The fit results for the NBO/FBO are shown in 
Figure 8 (see also Table 5). Below 52=2.1 the NBO (rep- 
resented by the filled circles in Figure 8) had a fairly con- 
stant frequency, with values between 6.3 Hz and 7.0 Hz. 
Its rms amplitude increased from 1.7% to 3.2%, and the 
FWHM was approximately ~2.1 or ~3.8. In the 52=2.0- 
2.1 selection both the NBO and FBO were present. This 
is likely an artifact of the 52 selection method, since a 
careful inspection of dynamical power spectra showed no 
evidence for simultaneous presence of both QPOs. It is 
interesting to note though, that in the 52=2.0-2.1 selec- 
tion the frequency of the FBO was 2.0±0.1 times that of 
the NBO, which could mean that the two NBO and FBO 
are harmonically related (however, see below). The FBO 
increased in frequency from 13.9 Hz (52 =2.04) to 23.1 Hz 
(52=2.65), while its FWHM increased from -5 Hz to -13 
Hz. The rms amplitude of the FBO initially continued the 
trend of the NBO rms amplitude; it increased from 2.8% 
(52=2.04) to 5.6% (52=2.23), but then decreased to 1.8% 
(52 =2.65). 



To study the transitions between the NBO and FBO 
more carefully, we inspected all dynamical power spectra 
of observations with 52 values around 2. Although no clear 
transitions were found, mainly due to the limited quality 
of the dynamical power spectra, wo did in some cases sec 
QPOs with intermediate frequencies (~ 10 Hz), suggest- 
ing that the frequency does not jump directly from ~7 
Hz to ~14 Hz. Fitting the power spectra in consecutive 
time intervals (longer than 16 s), rather than inspecting 
dynamical power spec:tra with the eye, will probably load 
to more conclusive results, but is beyond the scope of the 
current paper. The time scales on which the NBO/FBO 
frequency changed from —7 Hz to ~14 Hz and back were 
as short as a few tens of seconds. 

Wo also studied the behavior of the NBO/FBO during 
two long type I X-ray bursts. The first one started on 

1998 November 18 at 08:51:26 UTC, the second one on 

1999 Oct 10 at 09:10:47 UTC. Their exponential decay 
times were, respectively, 189±2 s and 144±2 s (Kuulkers 
et al. 2001). Both bursts occurred near the NB/FB ver- 
tex and in the power spectra of their respective observa- 
tions the NBO/FBO is clearly detected. No other QPOs 
were detected in the power spectra of these bursts. Fig- 
ure 9 shows the dynamical power spectra of both bursts, 
together with their 2-60 keV light curves. During the 
brightest part of the bursts the NBO/FBO seemed to dis- 
appear. Apparently, the burst flux was not modulated at 
the NBO/FBO frequencies with the same amplitude as the 
persistent flux. To quantify this, we determined upper lim- 
its on the NBO/FBO strength and compared those with 
the values outside the bursts. The results are shown in Ta- 
ble 6. Clearly, during the brightest part of the bursts the 
rms amplitude of the NBO/FBO was significantly weaker, 
not only as a fraction of the total fiux, but also in abso- 
lute terms (and hence as a fraction of the persistent flux, if 
the persistent emission is assumed to continue during the 
bursts). These measurements constitute the first determi- 
nation of the effect of X-ray bursts on the NBO/FBO. The 
fact that the bursts suppress the QPO can have important 
consequences for our understanding of its formation mech- 
anism (see Section 4.1.6). 

3.2.4. Noise components: LFN and VLFN 

LFN was detected between 52=— 0.6 and 52=5.0. As 

mentioned before, the appearance of the NBO around 
52=1.4, kept us from putting firm constraints on the LFN 
parameters between 52=1.4 and 52 = 1.6. The fit results 
for the LFN are shown in Figure 10 (see also Table 7), in 
two separate columns: one for the fits with a cut-off power 
law fit (52=0.0 1.4) and one for the fits with a Lorentzian 
(52<0.0 and 52=1.6-5.0). The strength of the LFN was 
in both cases defined as the integrated power spectral den- 
sity between 1 Hz and 100 Hz. We note that below 52=0.1 
the VLFN was not fitted separately from the LFN. An in- 
spection of 1/256-4096 Hz power spectra below 52=0.1 
showed that a weak power law component was present at 
frequencies below 0.1 Hz. This component was probably 
VLFN; its power in the 1-100 Hz range was much smaller 
than that of the LFN, so, although some VLFN power 
was absorbed by the LFN, this did not affect the LFN 
rms amplitudes significantly. The strength of the LFN 
changed considerably as a function of 52. It showed a 
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narrow peak between 5*2 =— 0.6 and Sz=0.2 of ~5% rms. 
Between 52=0.2 and 52 = 1.7 it gradually decreased from 
~4.5% rms to ^^2.5% rms. Another decrease was observed 
between 5^=2.7 and 5^=5.0, from ~3.8% rms to ~0.9% 
rms. The behavior between ^2=1. 7 and ^2=2. 7 was quite 
erratic, probably due to interactions with the fit functions 
of NBO and FBO. 

The centroid frequency of a Lorentzian and the cut-off 
frequency of a cut-off power law cannot be directly com- 
pared. Since we wanted to see how the typical frequency 
of the LFN evolved with Sz, we chose to plot Vmax, which 
is the maximum in a uP{v) plot and the frequency at 
which most of the power is concentrated per logarithmic 
frequency interval (Belloni, Psaltis, & van der Klis 2001, 
in prep.). For a Lorentzian v„iax is [vl + (FVFHM/2)2)i/2 
and for a cut-off power law Vmax is (1 — Oi)vcut (see Sec- 
tion 2 for analytical expressions for a Lorentzian and cut- 
off power law). As can be seen from Figure 10, between 
52=-0.6 and 52=1.7 

^max smoothly increased from ~1.2 
Hz to ~14.7 Hz. Above 52=1.7 the errors on Vmax were 
larger and the behavior was less clear. Between 52=1.7 
and 52=2.3 Vmax decreased to ^6.5 Hz, and above 52=2.3 
it increased again to ^15 Hz. We tested whether the 
change of the fit functions at 52=0.0 affected the values for 
^max and the other power spectral parameters, by swap- 
ping the fit functions used below and above 52=0.0; no 
significant changes were found. 

VLFN was detected over almost the whole 52 range. Al- 
though we only started fitting the VLFN separately from 
the LFN above 52=0.1 (in the 16 s power spectra), VLFN 
was present at frequencies below 0.1 Hz in the 256 s power 
spectra below 52=0.1. In the 16 s power spectra it was 
only significantly detected above 52=0.3. The fit results 
for the VLFN are shown in Figure 11 (see also Table 7). 
Between 52=0.3 and 52=1.4 the VLFN strength decreased 
from ^0.9% rms to ^0.4% rms. Above 52 = 1.4 its strength 
increased, to a peak value of ~1.3% rms at the NB/FB 
vertex. On the lower FB the strength decreased again, to 
a value of ^0.5% rms, and above 52=3.0 it increased to 
-1.6% rms at the top of the FB. The index of the VLFN 
slowly increased from —0.6 to —1.0 between 52=0.1 and 
52=2.4, with a small peak between 52=1.6 and 52=1.9, 
where the index had a value of ~1.5. Between 52=2.4 and 
52=3.2 the VLFN was much steeper. The indices were 
not well constrained and had values between 1.7 and 4.0. 
Above 52=3.2, where the error bars are much smaller, the 
index slowly decreased from 2.1 to 1.8. 

Our fits did not require an additional component to fit 
possible HFN; no such component was found. Possible 
explanations for this are discussed in Section 4.1.7. 

3.3. High Frequency QPOs 

Both the lower and upper kHz QPO were clearly de- 
tected; the lower kHz QPO between 52=0.5 and 52=1.5, 
the upper kHz QPO between 52=- 0.3 and 52 = 1.7. The 
results can be found in Table 8 and are shown in Figure 
12. The QPO frequencies showed a clear increase with 52, 
although both relations flattened at their low frequency 
ends. The FWHM of the lower kHz QPO was consistent 
with being constant at —100 Hz, whereas the rms ampli- 
tude showed a peak near 52=1.05 with a value of 3.6±0.3. 
The rms amplitude and FWHM of the upper kHz QPO 



both decreased with 52. 

Figure 13a shows the frequency difference of the two 
kHz QPOs as a function of the frequency of the upper kHz 
QPO. Some fits were made to the data; they are also shown 
in Figure 13a. The best fit to the frequency difference with 
a constant gave a value of 282±4 Hz. The ^ jd.o.j. for 
this fit was 18.6/9, which means that at a 97% confidence 
level the frequency difference was not constant. This is 
the first time that this is observed for GX 17+2. Fits 
with first and second order polynomials (taking into ac- 
count errors in both coordinates) resulted in, respectively, 
xVrf.o./.=15.0/8 and x^/d-o.f. =3.7/7. Although the lat- 
ter two fits show that with 99.76 % confidence (3.2(t) a de- 
crease towards higher frequency is not monotonic, it is not 
clear either whether the decrease towards lower frequencies 
is significant or not. This would be the first time that such 
a decrease towards lower frequencies is observed in any kHz 
QPO source with a non-constant frequency difference. To 
test this, two fits with a broken line (not shown) were per- 
formed, where in one case the slope of the low frequency 
part was fixed to zero. They resulted in /d.o.f. =6.28/7 
(slope fixed) and /d.o.f. =2.52/6 (slope free), which sug- 
gests that the frequency decrease towards lower frequency 
is only significant at a 97.6% confidence level (2.4cr). 

In Figure 13b we show additional fits to the data with 
theoretical curves for the radial epicyclic frequency vs. 
azimuthal frequency, which in the relativistic precession 
models are associated with the frequency difference and 
the upper kHz QPO frequency. The solid line is the best 
fit for circular orbtis (Stella & Vietri 1998; Markovic & 
Lamb 2001), which yields a xled °^ '^■^ (d.o.f = 9) for a 
mass of 2.07 Mq; the dashed line is the best fit for ec- 
centric orbits (Markovic & Lamb 2001), yielding a Xred 
1.14 (d.o.f = 8) for a mass of 1.86 Mq. Both curves were 
calculated for non-rotating stars; rotation of the neutron 
star gives only insignificant improvements. For compari- 
son we also plotted the data for Sco X-1 from van der Klis 
et al. (1997b). Clearly, both fits (especially the sohd line) 
do not fit the decrease of the frequency difference towards 
lower frequencies, which occurs in a frequency range for 
the upper kHz QPO that was not observed in Sco X-1, 
very well. 

The Q values (frequency/FWHM) of the two kHz QPOs 
were consistent with each other and with that of the HBO 
and its harmonic (see Figure 7); between 52 = 1.0 and 
52=1.5 the Q values of the HBO and kHz QPOh increased 
from —5 to —10. This unexpected finding may provide 
a key to understanding the formation mechanism of the 
three QPOs - this will be further discussed in Section 4.1.5. 

When comparing Figures 12 and 6 one can see that 
above 52=1.5, where the HBO frequency starts to de- 
crease, the upper kHz QPO frequency still increases. In 
Figure 14 we plot the frequencies of both QPOs against 
each other. For values of the upper kHz frequency lower 
than 1030 Hz the frequencies are well correlated, but above 
there is a clear anti-correlation. The solid line is the best 
power law fit to the points below 1000 Hz. The power 
law index is 2.08±0.07. This is the first time in any atoll 
or Z source showing kHz QPOs that an anti-correlation 
is observed between the frequencies of the low and high 
frequency QPOs (see Section 4.1.2 for a discussion). 
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4. DISCUSSION 

We performed a detailed study of the low and high fre- 
quency power spectral features of the neutron star low- 
mass X-ray binary and Z source GX 17+2. As was found 
in previous studies of GX 17-1-2 and other Z sources, the 
properties of most power spectral features correlated well 
with the position of the source along the Z track in the X- 
ray color-color diagram, suggesting that the spectral and 
variability changes are caused by one and the same param- 
eter. Some new results were found, the most interesting 
being the fact that the frequency separation of the kHz 
QPOs was not constant and probably not nionotonic (Fig. 
13), that their Q values were consistent with those of the 
HBO and its harmonic (Fig. 7), and that the frequency 
of the upper kHz QPO was anti-correlated with that of 
the HBO when the latter started to decline on the NB 
(Fig. 14). We also found a sub-harmonic of the HBO and 
showed that the absolute amplitude of the NBO/FBO is 
suppressed during type I X-ray bursts (Fig. 9). In the 
remainder of this section these and other results will be 
discussed in more detail. For a more detailed discussion, 
including a comparison of our results with those obtained 
from with other satellites (like EXOSAT and Ginga) and 
of other Z sources, we refer to Homan (2001). 

4.1. Timing behavior 

4.1.1. HBO and kHz QPO models 

In our discussion of the HBO and kHz QPO results, we 
will only consider models that try to explain the behav- 
ior of both types of QPOs. Currently these models are 
the relativistic precession model, the sonic point beat fre- 
quency model and the two-oscillator model. In the latter 
model (see Osherovich & Titarchuk 1999) an angle be- 
tween the magnetospheric equator and the disc plane is 
defined, which depends on the frequencies of the HBO 
and kHz QPOs and is supposed to be constant. How- 
ever, Jonker et al. (2001) show, by using our GX 17+2 
results and theirs of GX 5-1, that this is definitely not the 
case. For this reason, we will only consider the first two 
models. In the relativistic precession model (RPM, Stella 
& Vietri 1998, 1999; Morsink & Stella 1999; Stella et al. 
1999) three frequencies are predicted; they correspond to 
the fundamental frequencies of motion of a test particle in 
the vicinity of a compact object and are identified with the 
HBO and the two kHz QPOs. The lowest predicted fre- 
quency, which should explain the HBO, is a factor 2 lower 
than the low frequency QPOs in the atoll sources and the 
sub-harmonic of the HBO in the Z sources. Hence, a way 
to produce second harmonics is required. Note that no 
actual mechanism is proposed for the production of QPOs 
themselves. However, in a recent model by Psaltis & Nor- 
man (2001) the accretion disk acts as low band-pass filter 
with resonances close to the three frequencies predicted in 
the RPM. It is unclear to what extent the details of the 
RPM (i.e. the turnover in the relation between HBO and 
upper kHz frequency) carry over to the Psaltis & Norman 
(2001) description. We also note that recent investigations 
of the RPM by Markovic & Lamb (2001) have cast some 
doubt on the frequencies predicted by the RPM and re- 
lations between them. In the sonic point beat frequency 
model (SPBFM, Miller et al. 1998; Lamb & Miller 2001) 
the HBO is explained by the original beat frequency model 



of Alpar & Shaham (1985) and Lamb et al. (1985). The 
observed kHz frequencies are close to the orbital frequency 
at the sonic point of this disk and the beat of this frequency 
with neutron star spin frequency. The model requires the 
presence of a magnetosphere and a solid surface, and can 
therefore not explain the low and high frequency QPO in 
the black hole candidates. 

4.1.2. Horizontal Branch Oscillations 

The decrease of the HBO frequency on the NB in GX 
17-1-2 was first observed in RXTE data that were not used 
in our analysis (Wijnands et al. 1996). There were inter- 
esting features in the HBO properties which were repro- 
duced in that data and ours. In both cases the Sz depen- 
dence of the HBO frequency was not symmetric around 
its peak value; below, the frequency slowly flattened off to 
its maximum value, but above it showed a fast decrease 
by a few Hz that was followed by a more or less linear de- 
crease with Also, simultaneously with the frequency 
drop, the FWHM increased by a factor of more than two. 
Possible explanations for the HBO frequency decrease are 
discussed in Section 4.1.4. 

In all Z sources where an harmonic of the HBO has been 
observed, the frequency ratio of the harmonic and the fun- 
damental is generally less than two. This can be explained, 
at least in GX 17-1-2, by looking at the Sz dependence of 
the rms amplitudes of both components. In a certain Sz 
selection both frequencies tend to be dominated by those 
with the higher rms amplitudes, and hence by lower Sz 
and lower frequencies; this effect will be stronger for the 
harmonic than for the fundamental, since the Sz depen- 
dence of its rms amplitude is steeper (see Figure 6). As 
a result the frequency ratio will end up to be less than 
two. To test this hypothesis we calculated the weighted 
averages of two harmonically related linear functions, rep- 
resenting the HBO frequency-Szrelation between 5z=0.0 
and Sz—l-Q- As weight we used linear functions that 
were fit to the HBO rms-S'j; relation between 5*2=0.0 and 
5^=1.0. The ratio of the weighted frequency averages was 
1.93, which is very close to the observed value of 1.94. 

Peaked features with frequencies intermediate to those 
of the LFN and the HBO have been reported for GX 340-fO 
(Jonker et al. 2000b), Sco X-1 (van der Klis et al. 1997b; 
Wijnands & van der Klis 1999a), GX 5-1 (Kuulkers et al. 
1994; Jonker et al. 2001), and now also for GX 17+2. 
Based on the average frequency ratio of this sub-HBO 
component and the HBO in GX 17+2, 0.539+0.015, we 
suggest that, similar to what was suggested for GX 340+0, 
it may be the sub-harmonic of the HBO. A comparative 
study by Wijnands & van der Klis (1999a) of low frequency 
(HBO-like) QPOs and band limited noise components in 
the neutron star Z and atoll sources and black holes re- 
vealed strong correlations between the typical frequencies 
of those components. When plotting the frequency of the 
low frequency QPO versus the break frequency of the band 
limited noise they found that the branch traced out by the 
Z sources, although lying parallel to that of the atoll source 
and black holes, was slightly offset. It is interesting to see 
that this discrepancy between the Z sources on the one 
hand and the atoll sources and black holes on the other 
disappears (Jonker et al. 2000b), at least for GX 17+2, Sco 
X-1 and GX 340+0, when one uses the frequencies of the 
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sub-HBO component, instead of the HBO frequency (see 
also Wijnands & van der Klis 1999a). If the low frequency 
QPO in the atoll sources and the sub-HBO component of 
the three Z sources would indeed be the same QPO, this 
would mean that thc(sc) Z sources do not anymore follow 
the relation between the low frequency (HBO-like) QPO 
and the lower kHz QPO as defined by atoll and black hole 
systems (Psaltis et al. 1999). It is interesting to see that 
the three Z sources would then end up near or slightly 
above a second relation that is traced out by atoll sources 
like 4U 1728-34 and 4U 1608-52. At the moment, how- 
ever, it seems that the relation found by Wijnands & van 
der Klis (1999a) and the main relation of Psaltis et al. 
(1999) are mutually exclusive for the two types of sources. 
Note that Wijnands & van der Klis (1999a) found that if in 
their plot for Sco X-1 the frequencies of the sub-HBO com- 
ponent and the HBO were used, instead of the LFN and 
HBO frequencies, the discrepancy between Sco X-1 and 
the non-Z sources also disappeared. This reinterpretation 
has the advantage of preserving the consistency with the 
Psaltis et al. (1999) relation as well. 

The dominance of the even harmonics over the odd har- 
monics in the low frequency (HBO-like) QPOs, as sug- 
gested by the presence of a sub-harmonic of the HBO, is 
not only observed in the Z sources, but also a common fea- 
ture in BHCs such as GS 1124-68 (Belloni et al. 1997) and 
XTE J1550-564 (Homan et al. 2001). It suggests that a 
two-fold symmetry is present in the production mechanism 
of these low frequency QPOs. For the SPBFM one could 
think of an asymmetry in the magnetic field, which results 
in different areas for the polar caps (see e.g. Kuulkers et al. 
1997). In the RPM a two- fold symmetry is introduced by 
the two points at which the slightly inclined orbits of test 
particles cross the plane of the accretion disk. 

4.1.3. High frequency QPOs 

The kHz QPOs in GX 17+2 were discovered with RXTE 
(van der Klis et al. 1997a; Wijnands et al. 1997a). In our 
data the lower kHz QPO was found between 517 Hz and 
794 Hz (S'^=0.5-1.5) and the upper kHz QPO between 
618 Hz and 1087 Hz (S'2=-0.3-1.7), similar to what was 
reported by Wijnands et al. (1997a). Our findings of the 
kHz QPOs in GX 17-1-2 are compatible with the general 
properties of kHz in other Z sources (van der Klis 2000). 

An important new result is that the frequency differ- 
ence is no longer consistent (at a 97% confidence level) 
with being constant (see Figure 13). The frequency sep- 
aration reached a maximum value of 308 Hz, when the 
upper peak frequency was about 915 Hz, and fell off in 
both directions to values of 263 Hz and 239 Hz, when the 
frequencies of the upper peak were, respectively, 780 Hz 
and 1033 Hz. Note that the decrease in both directions, 
when fitted with a parabola, resulted in significantly (3.1 
cr) better fits than that with a constant. The decrease 
in the frequency separation at the low frequency side of 
the upper kHz QPO, although itself only significant at 
a 2.3(7 level, has not been seen before in any kHz QPO 
source. In the other LMXBs that show an non-constant 
frequency separation, Sco X-1 (van der Klis et al. 1997b); 
4U 1608-52 (Mendez et al. 1998a,b), 4U 1735-44 (Ford 
et al. 1998), and 4U 1728-34 (Mendez & van der Klis 
1999), the frequency separation was not measured at up- 



per peak frequencies as low as in GX 17-1-2. Although both 
the SPBFM and the RPM are capable of explaining the 
observed decrease in frequency difference towards higher 
upper peak frequencies, they have severe problems with a 
decrease towards lower frequencies. Such a decrease is not 
possible in the SPBFM, and, although it is in principle 
predicted by the RPM, that model only yields reasonale 
fits to our data when one allows for highly eccentric orbits 
(see Figure 13b). 

4.1.4. HBO-kHz QPO frequency relation 

The principal aim of our observing campaign in 1999 was 
to check on the relation between HBO frequency (whose 
increase was known to turn into a decrease towards higher 
Sz) and upper kHz QPO frequency. Clearly, a decrease, 
or even flattening, of the upper kHz QPO frequency, is 
not observed to accompany the decrease in HBO frequency 
(see Figures 6, 12, and 14). In a pure RPM, this can not be 
explained. However, if classical precession due to oblate- 
ness of the neutron star also plays a role, such an explana- 
tion is possible. For certain combinations of neutron star 
parameters, the prediction is that above a certain value of 
the upper kHz QPO frequency, the frequency of the HBO 
in fact starts to decrease (see Morsink & Stella 1999), due 
to the increasing importance of the classical precession. In 
GX 17+2 this turn-over of the HBO frequency occurs at 
a value of ~1 kHz for the upper kHz QPO with a peak 
frequency of ~60 Hz for the HBO. Judging from Figure 2 
in Morsink & Stella (1999), the relation between the HBO 
and upper kHz QPO frequencies in GX 17+2 would re- 
quire a high neutron star spin frequency (~500-700 Hz), 
a neutron star mass of more than 2 Mq, and a hard equa- 
tion of state (assuming the observed HBO frequency is 
four times the predicted frequency). In the RPM, well 
below the turn-over, the HBO frequency is expected to 
scale quadratically with the upper kHz QPO frequency. 
The best power law fit to the data below 1000 Hz yielded 
an index of 2.08+0.07 (see Figure 14), which is consistent 
with a quadratic relation. 

In the SPBFM a turnover in the HBO-upper kHz QPO 
relation is not predicted, but as the two frequencies do not 
directly depend on each other, the model is flexible enough 
to account for it; in case the turnover in HBO frequency 
can be explained by the two flow model put forward by 
Wijnands et al. (1996), the steady increase of the upper 
kHz QPO frequency could in principle be explained by as- 
suming that most of the radial flow falls back to the disk 
and rejoins the disk flow before the sonic point is reached, 
so that M at that location, and hence the upper kHz QPO 
frequency keep increasing. However, since the two radii at 
which the different QPOs are generated are very close, this 
seems a rather unlikely option. 

4.1.5. Q values of the HBO and kHz QPOs 

The Q values of the kHz QPOs are remarkably similar 
to each other, and to those of the HBO and its harmonic 
(see Figure 7). They also show an almost simultaneous in- 
crease above Sz=^, when the HBO harmonic is no longer 
detected. If the QPOs were artificially broadened by our 
selection method, averaging together peaks with different 
centroid frequencies, one would expect correlations of the 
FWHM (and Q value) with the width of the Sz selection 
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and with dv/dSz- Such correlations, if present at aU, are 
far too weak to explain the observed similarities. 

The FWHM of a QPO can be caused by several mech- 
anisms; intrinsic frequency variations, phase jumps, the 
simultaneous presence of several frequencies, and a finite 
lifetime of the signal will all cause a signal to appear as 
a broadened peak in a power spectrum. The data quality 
did not allow us to test the phase jump option. Lifetime 
broadening can probably not explain our results. In the 
SPBFM there is no reason why the lifetimes of the low 
and high frequency oscillations, which are determined by 
different processes, would lead to similar Q values. In the 
RPM, where all frequencies are generated by a single blob 
of matter, the lifetime for all variations should to first or- 
der be the same. However, in that case the Q values of 
the kHz QPOs should be much higher (by a factor 10-20) 
than those of the HBO, since more kHz cycles fit in a blob's 
lifetime. Of course, it is possible that the lifetime of the 
oscillations is not determined by the lifetime of the blob, 
but rather by a damping factor. Assuming that these fac- 
tors are not the same for the different QPOs, the similarity 
of the Q values could in principle be explained. However, 
there is no obvious reason why the (independent) damping 
factors would be fine tuned in a such a way that they lead 
to similar Q values. 

Explaining the Q values by broadening due to the pres- 
ence of several frequencies does not work either. Below 
the; turnover in the HBO-uppcr kHz QPO frequency rela- 
tion, the frequency of the HBO scales quadratically with 
the frequency of the upper kHz QPO. A given range of 
upper kHz frequencies would therefore lead to a Q value 
for the HBO that is half the Q value for the upper kHz 
QPO. This is not consistent with what we find. 

Frequency modulation could in principle explain the 
similar Q values, but only if the modulation is caused by 
a different parameter than the one that determines the 
frequency changes along the Z track. Frequency modula- 
tion only produces similar Q values if the QPO frequencies 
vary with the modulating parameter by the same factor, 
i.e. proportionally to each other. The observed quadratic 
relation between the HBO frequency and that of the up- 
per kHz QPO as well as the subsequent turnover, shows 
that this is not the case for the parameter that causes the 
motion along the Z track. However, assuming that the 
modulating parameter is not responsible for motion along 
the Z (see Section 4.2), we discuss the possibility of fre- 
quency modulation in a bit more detail. For broadening 
by frequency modulation to be observable in our power 
spectra, the time scale of these variations should be longer 
than that of the HBO. Additionally, the strength of those 
variations should anti-correlate with the Q values of the 
QPO. Assuming that the parameter that causes the fre- 
quency modulations also causes some changes in the count 
rate we could try to identify it in the power spectrum. 
The only reasonably strong components that we see in 
our power spectra and whose typical frequencies fulfill the 
above frequency requirement (v < 10 Hz) are the VLFN, 
the LFN and the NBO/FBO. The VLFN might be asso- 
ciated with the motion along the Z (see Section 4.1.7); 
disentangling that effect from a possible additional one on 
QPO frequencies is beyond the scope of this work. The 
strength of the LFN, which is of interest below 5^=1.5, 



anti-correlates with the Q values below S'z=1.5. Above 
Sz~l-5 the Q values decrease again to a value between 2 
and 5. In that case the broadening might be due to the 
appearance of the NBO/FBO and the increase in the LFN 
strength. It is interesting to note that in Sco X-1 it has al- 
ready been found that the kHz QPO properties are clearly 
modulated by the NBO (Yu et al. 2001). 

4.1.6. Normal and flaring branch oscillations 

For the first time we were able to study the proper- 
ties of the NBO/FBO during type I X-ray bursts. We 
found that the absolute rms amplitude significantly de- 
creased when the radiation from the neutron star surface 
increased. This is clearly different from the behavior of the 
~1 Hz QPOs in the dipping LMXBs 4U 1323-62 (Jonker 
et al. 1999), EXO 0748-676 (Homan et al. 1999), and 4U 
1746-37 (Jonker et al. 2000a), whose absolute rms ampli- 
tudes increased during type I X-ray bursts. It is thought 
that these QPOs are not related to the NBO/FBO in the 
Z sources. The behavior of the NBO/FBO in GX 17-1-2 
shows that NBO/FBO mechanism is very sensitive to the 
radiation field, thereby lending strong support for radi- 
ation feedback mechanisms (see below). In the frame- 
work of such mechanisms it is easy to understand that 
the NBO/FBO disappears during bursts: the NBO/FBO 
is only observed in a narrow Sz range, suggesting that it 
requires a delicate balance of certain parameters, among 
which the radiation field. The observations during the 
burst tails indicate that the NBO/FBO mechanism does 
not switch on instantly but rather becomes stronger in a 
more gradual way. Spectral fits during the bursts (Kuul- 
kers et al. 2001) suggest that the spectral properties of the 
accretion fiow itself are not affected much by the increase 
in radiation. 

The frequency of the NBO/FBO changed from being 
more or less constant below Sz=2 to being strongly cor- 
related with Sz, above Sz—2. In our data set we found 
the NBO between Sz=l-Q and 5^=2.1, although it may 
already have been present as a very broad feature as early 
as 5*^=1.4. Between 5^=2.0 and 5*^=2.1, the NBO (6.86 
Hz) was detected simultaneously with the FBO (13.9 Hz). 
Above Sz='2.l the FBO frequency increased, to ~23 Hz at 
5*^=2.65. Note that the simultaneous presence of the NBO 
and FBO between 5^=2.0 and Sz='2.1 is likely an artifact 
of the Sz selection method. Inspection of dynamical power 
spectra showed that they were never present at the same 
time. It is interesting though, that in that Sz selection 
their frequencies differed by a factor of 2, which is sugges- 
tive of mode switching of the NBO/FBO. However, there 
is evidence for QPOs with intermediate frequencies, which 
undermines that idea. The only way to see whether the 
frequencies really gradually transform into each other is 
to study the power spectral changes as a function of time. 
However, the quality of the GX 17-1-2 data is not sufficient 
to perform such a study. 

The presence of the NBO/FBO in the Z sources has of- 
ten been related to the fact that these sources accrete at 
near Eddington mass accretion rates. It is thought that 
at these high mass accretion rates a significant fraction of 
the accretion fiow is in the form of a thick and perhaps 
spherical fiow. Also, the effects of radiation pressure are 
thought to play an important role in this regime. The 
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model for the NBO by Fortner et al. (1989) is based on 
a radiation feedback mechanism in a spherical accretion 
flow. Another model for the NBO was proposed by Alpar 
et al. (1992), in which the NBO frequency is basically that 
of sound waves in a thickened accretion disk. This model 
does not explain how the NBO changes into the FBO. A 
major problem for both models could be the recently dis- 
covered NBO/FBO-like QPOs in the two atoll sources 4U 
1820-30 (Wijnands et al. 1999) and XTE J1806-246 (Wij- 
nands & van dcr Klis 1999b). Certainly in the first source 
the luminosity is believed to be only ~20-40 % of the Ed- 
dington luminosity, while near-Eddington luminosities are 
required in both models. 



4.1.7. Noise Components - LFN, VLFN 

The lack of HFN in our data, compared to, e.g., in the 
EXOSAT (Kuulkers et al. 1997) data can be explained 
by instrumental effects and different approaches to fitting 
the data. A study of EXOSAT d&ta, of several sources by 
Bcrger & van der Klis (1994) showed that a component 
with an rms amplitude of ~3% is always present in the 
frequency range where the HFN was often found in the 
EXOSAT data. When HFN was found at lower frequen- 
cies, it was, based on strength and frequency, most likely 
the component that we identified as the sub-harmonic of 
the HBO. 

Past measurements of the energy spectrum of the VLFN 
suggest that it is caused by motion of the source along 
the Z (van der Klis 1986, 1991; Lewin et al. 1992). The 
strength of the VLFN (Fig. 11) and the speed along the Z 
track {{Vz), see Figure 3) should therefore have a similar 
dependence on Sz- Our results are inconclusive; although 
some correlations are observed (the increase above iS'2=3, 
and the steep increase around Sz=2) there are also some 
differences. Most notable is the dip in the VLFN strength 
around 5*2=3, which is not observed as a decrease in (Vz). 
However, it is possible that due to the steepness of the 
VLFN around Sz=2.5 most of its power was outside the 
range in which we determined the strength. 

The LFN foflowed the behavior of the HBO, both in 
strength and frequency, until the latter disappeared {Sz = 
2.1). Between Sz='2.1 and Sz—^.S the behavior is rather 
complex and probably considerably influenced by the pres- 
ence of the NBO/FBO. Above 5^=3.0 the behavior was 
rather clean, with a decrease in strength and an increase 
in frequency. It is not clear, however, whether the LFN on 
the FB is the same as that on the HB and NB. Unfortu- 
nately no HBO was observed above Sz=2A to compare its 
behavior with. Van der Klis et al. (1997) suggested that 
in Sco X-1 the NBO/FBO emerged from the LFN. This 
seems not to be the case in GX 17-1-2, since the rms ampli- 
tude of the LFN already started to decrease well before the 
NBO appeared. Moreover, Vmax shows a gradual increase 
between Sz=-0.6 and Sz = l-4: to values well above the NBO 
frequency. A comparative study by Wijnands & van der 
Klis (1999a) strongly suggests that LFN has a similar ori- 
gin as the band-limited noise component in atoll sources, 
BHCs, and the millisecond X-ray pulsar SAX 1808.4-3658. 
They suggest that this noise component is produced in the 
accretion disk, and does not require the presence of mag- 
netosphere or solid surface. 



4.2. Broad-band Spectral Behavior 

In our observations GX 17-1-2 traced out the well known 
Z-like tracks in the CD and HID. In addition to motion 
along this track, we also observed secular motion of the 
Z track itself in the HID of epoch 3 (see Figure 1), as 
was already reported by Wijnands et al. (1997a) (see also 
Kuulkers et al. 1997, for a possible shift in the EXOSAT 
HID of GX 17+2). Below Sz = l the fluxes in the hard and 
soft energy bands behaved quite differently (Fig. 2); unlike 
the soft flux, the hard flux continued the trend that is ob- 
served above Sz~l- Hence, the HB/NB vertex in the CD 
is entirely due to the turn over of the soft spectral flux. 
It is not clear what causes the soft spectral flux to de- 
crease below Sz=l- Although it is partly an instrumental 
effect - the peak of the soft component slowly shifts out 
of the PCA effective energy range as it becomes softer- 
most likely intrinsic source changes occur at Sz=l- Most 
variability components increase in strength below Sz=l, 
suggesting that they are related to the hard spectral com- 
ponent. 

It is commonly bc;lievc!(l that the motion of the source 
along and its position on the Z track are determined by 
the mass accretion rate (M) through the inner accretion 
disk and onto the neutron star, increasing from the HB 
to the FB. Even though the count rates actually decrease 
on the NB and in some sources also on the FB (implying 
an anti-correlation between M and flux) there are several 
observational results that support this view: (1) the fre- 
quencies of the HBO and kHz QPOs, that in many models 
are directly related to M, gradually increase from HB to 
NB (and in Sco X-1 even to the FB, see van der Klis et al. 
1996). (2) The optical and UV fiux increase from HB to 
FB (Hasinger et al. 1990; Vrtilek et al. 1990; van Paradijs 
et al. 1990; Vrtilek et al. 1991; Hertz et al. 1992; Augusteijn 
et al. 1992); in X-ray binaries they are thought to be due 
to reprocessing of X-rays from the central source in the 
outer accretion disk. There are several reasons to believe 
that the optical and UV flux are betters tracers of the mass 
accretion rate than the X-ray flux (Hasinger et al. 1990). 
This would then solve the apparent anti-correlation be- 
tween M and the (X-ray) flux (Hasinger et al. 1990). (3) 
Motion of a source is always along the Z track; jumps be- 
tween branches are never observed, in accordance with the 
assumption that M varies continuously. 

It has been found that secular motion of the Z tracks in 
the CDs of the Cyg-like sources does not affect the relation 
between the rapid variability properties and source posi- 
tion along the track (see, e.g. Kuulkers et al. 1994; Jonker 
et al. 2000b, however see Wijnands & van der Khs 2001). 
In view of the M-driven picture mentioned above, these 
changes can therefore not be due to changes in M, but 
should have a different origin. It has been suggested that 
the Cyg-like sources are viewed at a higher inclination than 
the Sco-like sources (Kuulkers et al. 1994). The differences 
in shapes of the Z tracks and some of the variability prop- 
erties are also explained in this model (Kuulkers et al. 
1994; Kuulkers & van der Khs 1995). Psaltis et al. (1995) 
suggest that the difference among the Z sources might be 
explained by a difference in the magnetic field strength, 
with the Sco-like sources having a lower magnetic field. 

In the last few years results have been obtained that may 
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challenge the above described M-driven picture. These 
results are mainly due to the arrival of RXTE and Bep- 
poSAX, whose high quality data revealed new phenom- 
ena and allowed for a better comparison between different 
classes of X-ray binaries. In the following we discuss how 
some of our results of GX 17+2 conflict with the M-driven 
picture. 

The frequency of the HBO in GX17+2 starts to decrease 
when the source is half way through the NB. Although 
there are several explanations for this (see Section 4.1.2), 
it means that the first assumption given above (i.e. QPO 
frequencies increase with M) is not as solid as it seemed a 
few years ago. 

The theory of type I X-ray bursts predicts that their 
properties change and correlate with M (Fujimoto et al. 
1981; Fushiki & Lamb 1987). A study of the type I X-ray 
bursts in GX 17-1-2 by Kuulkers et al. (2001) showed that 
their properties did not correlate well with the position 
along the Z track at which they occurred, and therefore not 
well with the inferred M. Although other explanations are 
possible (see Kuulkers et al. 2001), this could mean that 
M does not determine the position of the source along the 
Z track. However, note that there does seem to be a corre- 
lation between burst occurrence and position along the Z 
track, which is usually explained by a changes in M onto 
the neutron star (along the Z track). 

Recently, di Salvo et al. (2000) reported the discovery 
with BeppoSAX of a hard tail in the energy spectrum of 
GX 17-1-2. It increased in strength from the HB/NB ver- 
tex (where only upper limits on its strength could be de- 
termined) onto the HB. They also showed that the total 
0.1-200 kcV flux increased monotonically from the NB/FB 
vertex (Fx = 1.52 x 10"® erg s^^) to the top of the HB 
{Fx = 1.84 X 10~® ergs~^). This is different from the 
behavior of the 2.9-20.1 keV count rate (see Figure 2a), 
which increased from NB/FB vertex to the HB/NB ver- 
tex, but decreases afterwards. As can be seen from Figure 
2b the behavior of the hard count rate confirms the find- 
ings of di Salvo et al. (2000). Also, there is no evidence 
for a change in its behavior near the HB/NB vertex. Un- 
fortunately no BeppoSAX data were taken on the FB to 
compare to our RXTE data. Hard power law tails have 
also been found in the Z sources GX 5-1 (Asai et al. 1994), 
Cyg X-2 (Prontera et al. 1998; di Salvo et al. 2001a), Sco 
X-1 (Strickman & Barret 2000; D'Amico et al. 2001) and 
GX 349+2 (di Salvo et al. 2001b). Except for Sco X-1, 
where no clear correlation with spectral state was found 
(D'Amico et al. 2001), in all these sources the strength 
of the hard power law tail seems to decrease in the same 
direction as was observed in GX 17+2, i.e. from HB to 
FB. 

The behavior of the 0.1-200 keV flux on the NB and HB 
suggest that M actually increases in the opposite direction 
from what is commonly assumed. However, since the flux 
changes are rather moderate (~20%), other options are 
that M does not change at all along the Z track, or at least 
that it shows no good correlation with the position along 
the Z. In that case motion along the Z track (i.e. spec- 
tral change) is caused by another parameter, for example 
the inner disk radius (van der Klis 2000). A similar pos- 
sibility was suggested by Homan et al. (2001) to explain 



the observed behavior of the black hole transient XTE 
J1550-564. They found that spectral transitions occurred 
at different levels of the inferred M. The spectral changes 
were; ascribed to an unknown parameter, which they sug- 
gested to be the inner disk radius (based on variability), 
the size of the Comptonizing region (based on spectral 
hardness), and/or the accretion flow geometry (based on 
radio brightness). Note that options like the latter one 
may lead to changes in the M onto the neutron star, while 
the M through the inner disk remains constant; in that 
way the occurrence of type X-rays bursts might become 
dependent on the position along the Z track. If the above 
also applies to the Z sources, the secular motion that is ob- 
served could perhaps be explained by changes in the M. 
The apparent lack of a correlation between the type I burst 
properties and the position along the Z track could then 
also be explained, since M would determine the properties 
of the bursts but no longer the spectral properties. 

It is interesting to see that GX 17+2 and the Z sources, 
show in fact several similarities with XTE J1550-564 and 
other black hole sources. The first was already briefly 
mentioned above, and concerns the independent behavior 
of the hard and soft flux on the HB. This independence 
of the hard and soft spectral components is commonly 
observed in black hole candidates (BHCs, see Tanaka & 
Lewin 1995). Moreover, there are several properties that 
correlate well with the spectral hardness in BHCs and also 
with the spectral hardness in Z sources, which increases 
from the NB to the HB (for the moment we will exclude 
the FB from the comparison; it will be discussed later). 
These properties are the radio flux, the optical flux, and to 
a certain degree also QPO frequencies. (1) The radio flux 
in BHCs is correlated with the spectral hardness (Fender 
et al. 1999; Fender 2001, and references therein). This is 
also true for GX 17+2, where the radio flux increases from 
the NB to the HB (Penninx et al. 1988). Sco X-1 (Hjellm- 
ing et al. 1990b) and Cyg X-2 (Hjellming et al. 1990a) show 
similar behavior, whereas GX 5-1 only showed one radio 
flare on the NB (Tan et al. 1992). The behavior of the 
radio flux suggests that as the spectral hardness increases, 
part of the flow onto the neutron star evolves into an out- 
flow. (2) Jain et al. (2001) found that in XTE J 1550-564 
the optical flux decreased during several small hard X- 
ray flares. Although the optical counterpart of GX 17+2 
(Callanan et al. 1999; Deutsch et al. 1999) has never been 
studied in much detail, this behavior is similar to the anti- 
correlation between the hard color and optical/UV flux on 
the NB and HB in the Z sources Cyg X-2 (Hasinger et al. 
1990; Vrtilek et al. 1990; van Paradijs et al. 1990) and 
Sco X-1 (Vrtilek et al. 1991; Hertz et al. 1992; Augusteijn 
et al. 1992). (3) In XTE J1550 564 the frequencies of the 
low and high frequency QPO anti-correlate with the hard 
flux (Sobczak et al. 2000; Homan et al. 2001). This is also 
found in GX 17+2 for the kHz QPOs and HBO (compare 
Figure 2b with Figures 6 and 12) and other Z sources. The 
only exception to this, together with the odd behavior of 
the hard tail in Sco X-1 (see above) is the correlation be- 
tween the HBO frequency and spectral hardness that we 
observed on the lower NB. In that respect it is interesting 
to mention that Rutlcdgc et al. (1999) report a switch from 
a correlation to an anti-correlation between the QPO fre- 
quency and spectral hardness in several BHCs, when the 
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latter increases beyond a certain value. In addition to the 
QPO frequencies also the band-limited noise behaves sim- 
ilarly. It becomes stronger with spectral hardness, both in 
BHCs and Z sources (van der Klis 1995a) 

It is tempting to compare the different black hole states 
(Tanaka & Lewin 1995; van der Klis 1995a) with those of 
the Z sources. Based on the above similarities, we suggest 
that the HE corresponds to the low state, the NB to the 
intermediate and/or very high state state, and the FB to 
the high state. However, there are still considerable differ- 
ences: the low state spectra are much harder than those 
of the HB, and also the noise in the low state is much 
stronger. Moreover, flaring is not observed in the black 
hole high state. However, in both the high state and FB 
the variability at high frequencies (>50 Hz) is extremely 
weak. It is not clear if and how the presence of a solid 
surface can account for these differences. 

5. SUMMARY AND CONCLUSIONS 

Our most important findings are the variable frequency 
difference of the kHz QPO, the non-monotonic relation 
between the HBO and the upper kHz QPO, the similar 
Q values of the HBO, its second harmonic, and the kHz 
QPOs and the suppression of the NBO amplitude during 
type I X-ray bursts. 

In the following we summarize what the implications 
of each of our new findings are for some of the proposed 
models. The models we considered for the HBO and kHz 
QPOs were the relativistic precession model (RPM), the 
'disk-filter' model, the sonic point beat frequency mode 
(SPBFM), and the two-oscillator model. The latter model 
can be discarded on the basis of our results (see Jonker 
et al. 2001). For the disk-filter model we assmned that 
it produces the same frequencies as the RPM and it will 
therefore not be discussed it separately. 

(1) HBO-kHz QPO frequencies: The RPM can only ex- 
plain a turnover of the HBO-kHz QPO frequency relation 
when one assumes that the observed HBO is four times the 
predicted frequency and then only for rather extreme neu- 
tron star parameters: a high neutron star spin frequency 
(^500 700 Hz), a neutron star mass of more than 2 Af©, 
and a hard equation of state. The SPBFM needs some 
extension; for example, decoupling of part of the accretion 
flow from the disk flow outside the magnetosphere (to ex- 
plain the HBO frequency decrease) and a subsequent re- 
coupling of this flow with the disk flow before the sonic 
point radius is reached (to explain the increase in the kHz 
QPO frequencies). How this should exactly work is not 
clear. 

(2) Variable kHz QPO separation: Although not signif- 
icant at a 3 sigma level, the kHz QPO separation seems to 



decrease both towards low and high values for the upper 
kHz QPO frequency. The RPM basically predicts such be- 
havior, but in its current form it is inconsistent with our 
results, as is apparent from Figure 13. The SPBFM can 
not explain the possible decrease of the frequency separa- 
tion at the low frequency side. 

(3) Similar Q values: In the SPBFM the low and high 
frequency QPOs are produced independently from each 
other and there is no a priori reason why their Q values 
should be similar. In the RPM, all the oscillations are 
produced by a single test particle. Lifetime broadening 
and the presence of multiple test particles both lead to 
different Q values for the low and high frequency QPO. 
Frequency modulation could work for both models, but 
only if the modulating parameter is not the same as the 
one that causes the motion along the Z track. 

It is clear that both the RPM and SPBFM in their cur- 
rent forms cannot satisfactorily explain all the observed 
phenomena. 

The observed behavior of the NBO during two type I X- 
ray bursts showed that this QPO is rather sensitive to the 
strength of the radiation field from the neutron star. This 
seems to confirm the basic idea behind the feedback mech- 
anism model that was proposed by Fortner et al. (1989), 
i.e. it requires a delicate balance between the inwards di- 
rected forces and the radiation field. 

Finally, we compared the behavior of GX 17-1-2 (and the 
other Z sourc;es) with that of (some) black hole LMXBs. 
Many similarities exist, particularly with regard to the re- 
lation between the hard spectral component on the one 
hand, and variability properties (QPO frequency, broad 
band noise strength) and radio flux on the other. Based 
on recent findings in black hole LMXBs we conclude that 
the mass accretion might not be the parameter that de- 
termines the position along the Z track after all. 
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Fig. 1. — Color-color diagrams (left column) and Hardness- Intensity diagrams (right column) for the epoch 3 (top) and epoch 4 data 
(bottom) . Each point represent a 16 s average. T?he splines and vertices (white circles) that were used for the Sz parameterization are shown, 
as are the typical error bars. See Table 3 for the energy bands used for the soft and hard color and intensity. 



15 




Fig. 2. — The count rate in three energy bands as a function of Sz (epoch 4 data only). The HB/NB vertex in the 2.9—20.1 keV band (a) 
is caused by the contribution from the low energies (b). This vertex does not show up at high energies (c). 
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17 




Frequency (Hz) 

Fig. 4. — Example of a fit to the 5z=0.4-0.5 power spectrum, showing the relative contribution of the individual components. The combined 
fit is represented by the thick solid line, the three Lorentzians/QPOs by the thin solid lines, the cut-off power law/LFN by the dashed line, 
and the power law/VLFN by the dotted line. The Poisson level was subtracted. The reduced x"^ of this fit was 1.03 for 204 degrees of freedom. 
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Fig. 5. — Power spectra (0.0625—256 Hz) for nine diflferent Sz selections. The Poisson level was subtracted for all power spectra. The most 
important power spectral features are indicated. 
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For reasons of comparison the frequency of the other two QPOs (gray) are also plotted in the frequency plot of each QPO. 
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Fig. 9. — The dynamical power spectra (top panels) and light curves in the 2-60 keV band (bottom panels) of the two type I X-ray bursts 
in which we studied the behavior of the NBO/FBO. The November 1998 burst is shown on the left (5 PCUs) and the October 1999, burst on 
the right (3 PCUs). The intervals in which we measured the NBO/FBO arc indicated by Roman numerals (see Table 6). The shades of gray 
in the dynamical power spectra represent the Lealiy power, with darker shades indicating higher powers. 
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FlK. 10. — LFN properties as function of Sz- The left column shows the results for the fits witfi a Lorentzian, the right column those for 
fits with a cut-off power law. In both cases the rms amplitude is the integrated power in the 1-100 Hz range. Umax is the frequency at which 
most of the power is concentrated (see text for the expressions for Vmax-) For reasons of comparison we also plotted the values of the other 
fit function (gray) in the panels for the rms amplitude and Vmax- 
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Fig. 11. — VLFN properties as function of Sz- The rms amplitude is the integrated power in the 0.1-1 Hz range. The arrows in the top 
panel represent upper limits. 
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Fig. 13. — The kHz QPO frequency diflfcrencc in GX 17+2 as function of the upper peak frequency. In (a) three fits to the data are shown 
(gray): a constant (solid line), a straight line (dashed line), and a parabola (dotted line). In (b) the two best fits to the data are shown, for 
theoretical relations between the radial epicyclic frequency and aziniuthal frequency for circular orbtis (solid line) and eccentric orbits (dashed 
line). The fits were kindly provided by Draza Markovic. Note that the dashed line requires highly eccentric orbits with apastron to periastron 
ratios of ~2-3. For comparison we also plotted the data for Sco X-1 (gray) from van der Klis et al. (1997b). For additional details see text 
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Fig. 14. — HBO frequency as function upper kHz QPO frequency. The solid line is the best power law fit to the data below 1000 Hz for 
the upper kHz frequency. The power law index is 2.08±0.07. A clear deviation from this relation can be seen for values above 1000 Hz. 
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Begin Time (UTC) 


End Time (UTC) 


Total (ks) 


Modes" 


Gain Epoch 


1997-02-02 19:13 


1997-02-27 03:34 


58.7 


3 or 4,5,6,9 


3 


1997-04-01 19:13 


1997-04-04 23:26 


34.6 


4,5,6,9 


3 


1997-07-27 02:13 


1997-07-28 00:33 


42.9 


4,5,6,9 


3 


1998-08-07 06:40 


1998-08-08 23:40 


71.0 


4,5,8 


3 


1998-11-18 06:42 


1998-11-20 13:31 


86.0 


4,5(,7),8 


3 


1999-10-03 02:43 


1999-10-12 07:05 


297.6 


10,11,12 


4 


2000-03-31 12:15 


2000-03-31 16:31 


6.9 


10,11,12 


4 



" Modes in addition to the Standard 1 and Standard 2 modes. See Table 2 for modes. 



Table 1 

a log of all rxte/pca observations used in this paper. mode 7 was not always active during the november 

1998 observations. note that none of the observations represents an uninterrupted interval. each is a 
collection of observations that were done around the same time. these observations were separated in time 
from each other for various reasons, such as earth occultations, passages of the south atlantic anomaly, or 

observations of other sources. 
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Mode 


Name 


Time Resolution (s) 


Energy range (keV) 


Energy channels 


1 


Standard 1 


2-3 


2 60 


1 


2 


Standard 2 


2^ 


2-60 


129 


3 


E_8us_8A_0_ls 


2-17 


2-60 


8 


4 


SB_125us_0_13_ls 


2-13 


2-5.1 


1 


5 


SB_125us.l4.17_ls 


2-13 


5.1-6.6 


1 


6 


SB_125us.l8.23_ls 


2-13 


6.6-8.7 


1 


7 


SB_125us.l8.249_ls 


2-13 


6.6-60 


1 


8 


E_16us_64M.18_ls 


2-16 


6.6-60 


64 


9 


SB_125us.24.249_ls 


2-13 


8.7-60 


1 


10 


SB_125us.0.13.1s 


2-13 


2-5.8 


1 


11 


SB_125us.l4.17_ls 


2-13 


5.8-7.5 


1 


12 


E_16us_64M.18_ls 


2-16 


7.5-60 


64 



Table 2 

Names and settings of the data modes that were used in our analysis. The lower and upper energy boundaries 
OF the PCA energy sensitivity range are given as 2 and 60 keV, although they changed between (and also 

slightly during) the different gain epochs. 
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Gain Epoch 


Soft Color 


Hard Color 


Intensity 


3 
4 


(Channels) (keV) 
10-16/5-9 4.8-7.3/3.0-4.8 
8-13/4-7 4.6-7.1/2.9-4.6 


(Channels) (keV) 
26-50/17-25 10.5-19.7/7.3-10.5 
22-42/14-21 10.5-19.6/7.1-10.5 


(Channels) (keV) 
5-50 3.0-19.7 
4-42 2.9-19.6 



Table 3 

Channel and energy boundaries of the soft and hard colors and the intensity used for the spectral analysis. 
The channel numbers refer to the Standard 2 mode channels (1-129). 
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Component 


S'x-rangc 


LFN 


-0.6 


-5.0 


HBO (fund.) 


-0.6 


-2.1 


HBO (2nd harm.) 


-0.6 


-1.0 


HBO (sub-harm.) 


-0.6 


-1.3 


Upper kHz QPO 


-0.3 


-1.7 


VLFN 


0.3- 


5.5 


Lower kHz QPO 


0.5- 


1.5 


NBO 


1.6 


2.1 


FBO 


2.0- 


2.7 



Table 4 

The nine different components in the combined epoch 3/epoch 4 power spectra, and the S^-ranges in which 
they were detected. the components are listed in order of appearance. 



HBO Fundamental 



HBO 2"'' Harmonic 



HBO Sub-harmonic 



NBO 



Rms" (%) 



I Frequency (Hz) Rms" (%) 



Frequency (Hz) Rms'' (%) 1 



I Frequency (Hz)Rms'' (%) FWHM (Hz) Frequency (Hz) 



— 0.43±0.05 


4.3±0.7 




21.3±0.8 


4.5±0.7 


9^3 


41.4±1.4 


6.4±1.0 


rtfr+lS 

27+1° 


(fixed) 








-0.32±0.04 


3.9±0.6 




21.4±0.4 


4.9±0.3 


11.7±1.7 


41.4±0.5 


6.9±0.4 


19±6 


(fixed) 








-0.23±0.03 


2.9±0.3 


4.1±0.9 


22.2±0.2 


3.9±0.2 


7.8±1.1 


43.9±0.3 


7.4±0.3 


44±6 


(fixed) 








— U. iuitU.Uo 


O.OltU.O 


? Q-i-1 n 

1 .Olt i.U 


ZZ.OZlZKJ.O 


t- 14+0.19 


1 c; Q_i_i A 




D. 1 itU..^ 




u [jixeaj 








-0.04±0.03 


4.02±0.20 


8.7±0.9 


24.3±0.2 


4.79±0.12 


1.5. Oil. 


47.3±0.3 


6.41±0.17 


19 + 3 


(fixed) 








0.05±0.03 


3.12±0.14 


5.4±0.5 


26.39±0.13 


4.06±0.14 


12.7±0.9 


50.7±0.2 


7.0±0.9 


29±3 


0±4 








U. 1 D ^ U. IJO 


Q 9 1 _i_n IS 


6 1 ih 6 


97 zl9-l-n 1 zL 


s9+n 1 


12 4i0 9 


53 7ibO 3 


D.D_0.9 

5.7it0.8 


31it4 


9 +3 








{).Zo±u.()o 


3. 54ih0. 17 


8.8±0.7 


30.20±0.18 


0. / o±U. iz 


16.4i 1 .2 


( .y±u.4 


25ib3 


8it3 








0.35±0.04 


2.81±0.11 


6.8±0.5 


33.33±0.15 


2.86±0.12 


l.'5.0±1.3 


63.6±0.4 




32^3 


n + 2 

— 8 








n /I fi-t-n ri/i 


2 95ih0 1 1 


. m u . 


3.^ 3'5i0 1.5 


1 no 


13 9 i 1 2 


Do. .1 


4.1±0.4 




-1-" 1,2 








0.56±0.03 


2.75±0.09 


8.3±0.5 


37.91±0.14 


2.41±0.08 


13.6±1.1 


73.4 ±0.3 


26±3 


15. Oil. 5 








0.65±0.03 


2.57±0.10 


7.8±0.6 


39..56±0.14 


2.42±0.09 


16.8±1.5 


77.9±0.5 


3.9±0.2 


28±3 


18. Oil. 








0.75±0.03 


2.41±0.13 


9.1±0.9 


42.9±0.2 


1.92±0.11 


14. ,^±1. 9 


83.3±0.7 


3.6±0.3 


28±4 


20.6±1.0 








0.85±0.03 


2.2±0.2 


9.8±1.6 


45.3±0.4 


1.62±0.1.5 


13±3 


89.4±0.9 


3.2±0.3 


29±5 


24.9±1.4 








0.95±0.03 


2.3±0.3 


13. ±3 


48.9±0.7 


1.8±0.2 


22I5 


93±2 


2.5±0.4 


21±7 


25.2±1.3 








1.06±0.03 


1.52±0.17 


6.0±1.3 


.■33.9±0.4 


1.09±0.15 


15 (fixed) 


107 (fixed) 


'^-'O_0.19 


31±7 


32.2±1.6 








1.15±0.03 


1.40±0.12 


5.1±1.0 


.■36.7±0.3 


<0.7 


15 (fixed) 


113 (fijced) 


9 r+0.3 


32±6 


35.5±1.6 








1.25±0.03 


1.42±0.12 


r 1+0.8 

°-l-0.3 

4.0±0.6 


58.5±0.3 








n Q+0.5 
^■*'-0.3 

< 2.5 


°' -12 


41±6 








1.35±0.03 


1.63±0.11 


60.11±0.17 








21±14 


43±3 








1.45±0.04 


1.40±0.12 


5.2±1.2 


60.3±0.3 




















1.55±0.03 


1.8±0.2 


14±3 


54.6±1.0 




















1.65±0.04 


1.72±0.14 


12±3 














1.7±0.2 


3.8±0.7 


6.33±0.16 


1.76±0.04 


1.99±0.12 


17±3 


53.6±1.0 














2.07±0.11 


2.08±0.20 


6.82±0.04 


1.86±0.04 


1.90±0.12 


17±3 


52.6±1.1 














2.89±0.08 


2.13±0.11 


6.91±0.03 


1.95±0.04 


1.80±0.16 


19±4 


50.1±1.2 














2.94±0.12 


2.11±0.13 


6.96±0.02 


2.04±0.04 


1.2±0.2 


llt^ 


48.2±2 














3.2±0.2 
2.8±0.8 


3.7±0.4 


6.86±0.10 


2.13±0.05 


<1.2 


9 (fixed) 


46 (fixed) 














4.93±0.11 


5.1±0.2 


14.78±0.08 


2.23±0.06 
















5.60±0.07 


6.4±0.18 


15.44±0.06 


2.34±0.05 




















5.28±0.11 


8.1±0.3 


16.75±0.12 


2.45±0.05 




















4.37±0.17 


10.3±0.7 


18.2±0.2 


2.54±0.04 




















3.36±0.20 


13.6±1.6 


19.7±0.6 


2.65±0.04 




















1.8±0.4 


8±4 


23.1±0.9 


2.75±0.04 




















< 1.8 


10 (fixed) 


25 (fixed) 



" Integrated between — oo Hz and +00 Hz. 

^ Integrated between Hz and +00 Hz, because of low Q values. 



Table 5 

Fit results for the low frequency QPOs. 



Int(Tval 




Count Rate 


Fractional rms 


Absolute rms 






(counts s^'-j 


(%) 


(counts s^"'") 


Nov 1998 - 


I 


2811 


2.58±0.19 


72.5±5.3 


Nov 1998 - 


II 


2648 


3.2±0.3 


85±8 


Nov 1998 - 


III 


4568 


<0.71 


<32.4 


Nov 1998 - 


IV 


3232 


<2.1 


<68 


Nov 1998 - 


V 


2833 


3.3±0.2 


93±6 


Oct 1999 - 


I 


1742 


5.56±.17 


96.9±3.0 


Oct 1999 - 


II 


2616 


<1.57 


<41.1 


Oct 1999 - 


III 


1847 


4.3±0.4 


79±7 



Table 6 

The properties of the NBO/FBO during the two bursts shown in Figure 9. The intervals given in first coluimn 

CAN also be found IN THAT FIGURE. THE UPPER LIIVIITS IN THE LAST TWO COLUMNS ARE 95% CONFIDENCE. THE COUNT 
RATES ARE IN THE 5.1-60 KEV (NOVEMBER 1998, 5 PCUs) AND 5.8-60 KEV (OCTOBER 1999, 3 PCUs) BANDS. 
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VLFN 


Ll'N (Cut-ott power 


law) 




Ll'N (Lorontz 


an) 




Rms" (%) 


Index 


Rms" (%) 


Index 


OUt-OII l^rlzj 


Kms { /q) 


r W rliVi jlz J 


Frequency (Hz) 
















Tf^^TTn^ — 

z . o i U . 


'1 zl+d 7 


i \ U+ll 4 

u.ymu.o 


-0.32±0.04 














3.7±0.9 


3.1±1.0 


1.05±0.19 


-0.23±0.03 














5.0±0.3 


4.0±0.4 


1.15±0.13 


— 40.15±0.03 














4.5±0.5 


3.9±0.6 


1.19±0.09 


-0.04±0.03 








■ ■ ■ 






4.2±0.4 


3.4±0.4 


1.30±0.06 


0.05±0.03 








4.1itU.3 


n 'vc _i_n 1 1 
— U. /DihU.ii 


1.43±0.15 












Q 


U.D [JTiXGdJ 


4l:,oOiU. 10 


n sri-l-n 1 1 










0.25±0.03 




<0.7 


0.6 (fixed) 


4 6i0 3 




1 93±0 19 








0.35±0.04 





86±0.07 


0.69±0.04 


, ,0 + 0-14 

4.46^0.24 


— 0.86±0.07 


1.85+"!," 

— 0.08 








U.4b±U.U4 







0.61±0.03 


4.44+0-13 


-o.86tj;°;^ 


9 ni +0.06 








0.56±0.03 





DO±U.ub 


0.66±0.0,'j 


4.66±0!ll 


— 0.71±0.05 


2.54±0.14 








0.65±0.03 





65±0.06 


0-85_o.o9 
U.o7±U.Uo 


4. ( 4±U.U / 


— U.u ( ±U.U4 


2.94±0.12 








0.75±0.03 


u 




4.oo±U.U7 


— U.7U±U.Uo 


3.24±0.14 








0.85±0.03 





61±0.08 


0.47±0.05 


4.53±0.08 


-0.61±0.06 


4.0±0.2 








0.95±0.03 





67±0.07 


0.64+00^ 

— U.Oo 

0.4 (fixed) 


4.26±0.12 


-0.58±0.08 


4.8±0.4 








1.06±0.03 




<0.5 


4.04±0.13 


-0.33±0.05 


6.9±0.6 








1.15±0.03 





60±0.05 


0.58±0.04 


3.52±0.07 


— 0.07 


6.8±0.4 








1.25±0.03 




<0.5 


0.5 (fixed) 


3.55±0.12 


-0.156_o 021 


— O.o 








1.35±0.03 





43±0.09 


0.60±0.15 


3.24±0.11 


-0.29±0.10 


9.2±1.2 








1.45±0.04 





63±0.04 


„ „„ + 0.10 

0.54±0.03 














1.55±0.03 





71±0.04 














1.65±0.04 





62±0.05 


1.29±0.19 








2.2±0.3 


2i±5 


lb±2 


1.76±0.04 





62±0.04 


1.6 ±0.2 








2.55±0.13 


15.4±1.8 


9.2±0.7 


1.86±0.04 





70±0.03 


1.48±0.14 








2.63±0.12 


13.5±1.3 


9.9±0.6 


1.95±0.04 





82±0.03 


0.92±0.10 








3.42±0.17 


9.4±0.3 


9.6±0.3 


2.04±0.04 


1 


34±0.03 


0.77±0.03 








2.9±0.9 


5.3±1.3 


11.0±0.5 


2.13±0.05 


1 


03±0.05 


l.OlibO.ll 








3.6±0.2 


7.1±0.8 


7.1±0.20 


2.23±0.06 





85±0.05 


86+° " 








2.50±0.17 


6.0±0.9 


5.8±0.2 


2.34±0.05 


1 


02±0.08 


1 ,,+0.23 












5.6±0.7 


2.45±0.05 





52±0.05 










^■"-0.3 


12tl 


3 9+° '' 


2.54±0.04 





67±0.06 


2.5±0.4 








2.7±0.2 


9. Oil. 6 


4.8±0.4 


2.65±0.04 





66±0.08 


1.7±0.4 








3.5±0.3 


16±3 


5.8±0.7 


2.75±0.04 





59±0.07 


9 O + 0.9 








3.82±0.17 


19±3 


5.9±0.4 


2.85±0.04 





46±0.06 


q 9+1.5 
■^■^-0.8 

2.7±0.5 








3.48±0.14 


14.7±1.6 


5.3±0.4 


2.95±0.04 





55±0.05 








3.66±0.11 


15.9±1.3 


6.0±0.4 


3.10±0.06 





64±0.03 


2.6±0.3 








3.42±0.06 


14.1±0.7 


6.1±0.2 


3.29±0.07 





80±0.03 


2.12±0.18 








3.23±0.09 


16.3±1.1 


6.3 ±0.4 


3.50±0.07 





85±0.03 


2.14±0.16 








2.84±0.08 


17.4±1.3 


7.0±0.5 


3.70±0.07 


1 


00±0.03 


1.83±0.10 








2.23±0.09 


15.7±1.7 


9.5±0.5 


3.89±0.06 


1 


03±0.03 


1.90±0.10 








1.75±0.09 


12.6±1.7 


10.7±0.7 


4.12±0.08 


1 


03±0.03 


1.90±0.10 








1.49t"„;- 




7.01-", 


4.37±0.08 


1 


22±0.03 


1.79±0.08 








i-4otS:r6 


23+^° 


9tl 


4.71±0.15 


1 


65±0.03 


1.81±0.06 








0.94± 0.12 


11 ±3 


8 (fixed) 


5.13±0.10 


1 


61±0.04 


1.80±0.08 








< 0.6 


11 (fixed) 


8 (fixed) 



" Integrated between 0.1 Hz and 1 Hz 
^ Integrated between 1 Hz and 100 Hz 



Table 7 

Fit results for the noise components at low frequencies. 







Lower kHz QPO 




Upper kHz QPO 






Rms (%) 


FWHM (Hz) 


Frequency (Hz) 


Rms (%) 


FWHM (Hz) 


Frequency (Hz) 


Difference (Hz) 


-0.32±0.04 








<7 








-0.23±0.03 








4.4±0.6 




637±9 




-0.15±0.03 








5 4+1-0 


192115 


618±12 




-0.04±0.03 








5.7±0.4 


174±25 


650±7 




0.05±0.03 








5.5±0.4 


189±29 


663±8 




0.16±0.03 








5.4±0.4 


178±26 


678±7 




0.25±0.03 








5.8±0.4 


205±24 


689±6 




0.35±0.04 








4.8±0.3 


175±21 


722±6 




0.46±0.04 


<2.0 


130 (fixed) 


475 (fixed) 


4.6±0.2 


164±16 


759±5 




0.56±0.03 


2.4±0.4 


145±55 


517±18 


4.4±0.2 


167±18 


780±5 


263±18 


0.65±0.03 


9 7 + 0-4 
^- ' -0-3 
1.6±0.4 


137lf„ 




4.46±0.20 


160±15 


804±4 


284±11 


0.75±0.03 


85±59 


530±20 


4.3±0.2 


181±21 


823±7 


293±20 


0.85±0.03 


2.4±0.3 


87±25 


555±8 


4.0±0.3 


154±20 


849 ±6 


294±10 


0.95±0.03 


3.1±0.3 


104±25 


607±8 


3.5±0.4 


156±33 


915±11 


308±14 


1.06±0.03 


3.6±0.3 


144±24 


655±8 


3.5±0.3 


149±30 


950±8 


294±11 


1.15±0.03 


2.9±0.2 


87±15 


682±5 


3.2±0.4 


200±49 


953±14 


271±15 


1.25±0.03 


, n+O-3 

^■^-0.2 


112tf 


714±7 


^■^-0.3 




973±11 


258±13 


1.35±0.03 


2.0±0.3 


67±24 


752±8 


1.8±0.4 




1011±24 


259±25 


1.45±0.04 


2.0±0.3 


98±41 


794±14 


1.7±0.3 


61±28 


1033±9 


239±17 


1.55±0.03 


<2.3 


80±59 


830±19 


1.8±0.3 


76±32 


1063±13 




1.65±0.04 








l-°-0.3 


83lf, 


1087±15 




1.76±0.04 








<1.5 


90 (fixed) 


1085 (fijced) 





Table 8 

Fit results for the high frequency QPOs. 
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